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ABSTRACT 
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The Microvariability and Oscillations of Stars (MOST) satellite observed the 
^ ■ B supergiant HD 163899 (B2 Ib/II) for 37 days as a guide star and detected 48 

Q\ • frequencies < 2.8 c d -1 with amplitudes of a few milli-magnitudes (mmag) and 

less. The frequency range embraces g- and p-mode pulsations. It was gener- 
^ ■ ally thought that no g-modes are excited in less luminous B supergiants because 

^ ■ strong radiative damping is expected in the core. Our theoretical models, how- 

ever, show that such g-modes are excited in massive post-main-sequence stars, in 
accordance with these observations. The nonradial pulsations excited in models 
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between 2OM at logT eff m 4.41 and 15M© at logT cff pa 4.36 are roughly consis- 
tent with the observed frequency range. Excitation by the Fe-bump in opacity 
is possible because g-modes can be partially reflected at a convective zone as- 
sociated with the hydrogen-burning shell, which significantly reduces radiative 
damping in the core. The MOST light curve of HD 163899 shows that such a 
reflection of g-modes actually occurs, and reveals the existence of a previously 
unrecognized type of variable, slowly pulsating B supergiants (SPBsg) distinct 
from a Cyg variables. Such g-modes have great potential for asteroseismology. 

Subject headings: stars: individual (HD 163899) - stars: early-type - stars: 
oscillations - supergiants 

1. INTRODUCTION 

The B star region in the HR diagram contains a wide range of variables. The most 
luminous (luminosity classes of O-Ia; > 1O 5 L ) are known as a Cyg variables and show semi- 
regular light and radial- velocity variations on timescales of 1 - 6 weeks (e.g., van Leeuwen 
et al. 1998). Those variations are thought to be a combination of pulsation and rotational 
modulation of winds (e.g., van Genderen 1991; Kaufer et al. 1996, 1997). Pulsations in such 
very luminous stars can be excited mainly by the Fe-bump of opacity (Rogers & Iglesias 1992) 
at a temperature of about 2 x 10 5 K. The excitation mechanism is, however, more affected by 
the strange-mode character in more luminous stars (Kiriakidis et al. 1993; Glatzel & Mehren 
1996; Dorfi & Gautschy 2000). The target star discussed in this paper, HD163899 (B2Ib/II), 
is a supergiant that is less luminous and less affected by winds than the a Cyg variables. 

There are also many less luminous B-type variable stars whose light and radial-velocity 
variations are caused by radial and nonradial pulsations. These pulsations are excited by 
the classical K-mechanism at the Fe-bump of opacity. The type of mode excited depends 
mainly on the effective temperature of the stars; radial and nonradial p-mode pulsations are 
excited in hotter stars and nonradial g-mode pulsations are excited in relatively cooler stars 
(see e.g., Pamyatnykh 1999). Most of the stars are main-sequence stars. Hotter p-mode 
pulsators {(3 Cephei stars) are more luminous or more massive (> 8M Q ) than the relatively 
cooler (B3 - B8) g-mode pulsators, slowly pulsating B (SPB) stars (Waelkens 1991). Mainly 
based on line-profile variations, g-modes are known to be excited also in early-type Be stars 
(e.g., Rivinius et al. 2003). 

The (3 Cephei stars have been known for more than a century, but the excitation mech- 
anism was identified only in the early 1990s (Kiriakidis et al. 1992; Moskalik & Dziembowski 
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1992) after the emergence of OPAL opacity tables (Rogers & Iglesias 1992). The obser- 
vational properties of the (3 Cephei stars are thoroughly discussed by Stankov & Handler 
(2005). 

The SPB stars are nearly perfectly confined to the main- sequence band (e.g., De Cat et 
al. 2004), in accordance with the theoretical prediction that the radiative core in a post-main- 
sequence star strongly damps g-modes where they have very short wavelengths (Gautschy 
& Saio 1993; Dziembowski et al. 1993). In this paper, however, we show that reflection of 
pulsation modes at the convective zone associated with the hydrogen-burning shell quenches 
the radiative damping in the core, and hence the Fe-bump can excite some g-modes in a 
post-main-sequence massive star with a radiative core. 

The Microv (inability & Oscillations of S Tars (MOST) photometric satellite was launched 
on June 2003 and is fully described by Walker et al. (2003). The first scientific results were 
published by Matthews et al. (2004). 

MOST observations have significantly clarified the properties of pulsating B stars. 
MOST detected clear light variations in the rapidly rotating 09Ve star ( Oph (Walker 
et al. 2005a), which revealed the presence of (3 Cephei- type pulsations in this line-profile 
variable. MOST found that the mono-periodic (3 Cephei star S Ceti is actually multiperiodic 
with three additional frequencies (Aerts et al. 2006a). MOST discovered a new SPB star, 
HD 163830 (B5II/III), near the edge of the main-sequence band, detecting 20 g-mode fre- 
quencies (Aerts et al. 2006b). MOST detected g-modes and possibly r- modes in the rapidly 
rotating Be star HD 163868 (Walker et al. 2005b). The latter observations showed, for the 
first time, that many g-modes can be excited in a rapidly rotating Be star, and provided 
a natural interpretation of A Eri-type variables as SPBe stars. This finding suggests that 
nonradial pulsations might play a crucial role in Be phenomena. 

In this paper, we report another important discovery by the MOST satellite - g-modes 
in the B supergiant HD 163899. 

2. THE MOST OBSERVATIONS 

HD 163899 (V=8.3, B2 Ib/II) was observed as one of 20 guide stars for the photometry 
of WR 103 (HD 164270) by the MOST satellite. The observing conditions and the data 
reduction method were the same as for the SPBe star HD 163868 described in Walker et al. 
(2005b). Observations were made from 2005 June 14 to July 21 for a total of 36.6 days with 
> 137, 000 values recorded for HD 163899. The data are binned in 2 minute intervals. Since 
the WR 103 field is outside the MOST Continuous Viewing Zone, the duty cycle was limited 
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to about 50% of the 101 min orbit. 

Figure 1 displays the full light curve of HD 163899 and for another guide star, HD 
164388 (A2 IV, V=8.0), which appears to be constant (standard deviation ~ 1.6 mmag). 
An expanded two day portion of the MOST light curve is included to better display the 
higher frequency variations. The light curve shows many periodicities ranging from about 
a day to a few days, indicating that many g-mode pulsations are excited in HD 163899. 
Complete light curves for both stars can be downloaded from the MOST Public Archive at 
www.astro.ubc.ca/MOST. 

2.1. The Frequency analysis 

The frequency analysis was carried out using CAPER, a collection of Fortran driver 
routines which use a Discrete Fourier Transform (DFT) as a frequency and amplitude iden- 
tification tool and non-linear least squares (NLLS) software (Numerical Recipes; Press et al. 
1986) to refine the identified parameters. The fitting function is a series of sinusoids (see 
Walker et al. 2005b). 

Periodicities were identified and refined as follows: 

a) A frequency/amplitude was identified from the DFT. The amplitude and phase were 
refined via NLLS. The fit function was subtracted from the data and a new set of parameters 
was identified from the residuals. The new amplitude and phase was refitted along with all 
other previously identified amplitudes and phases. The new function was then subtracted 
from the original time series and new parameters identified. This process is repeated until 
the next identified frequency has a S/N ~ 3.6 (> 3 sigma detection) (Kuschnig et al. 1997). 

Noise is defined as the mean of the amplitude spectrum ±3 c d _1 centered on the 
identified frequency. The amplitudes are 3 sigma clipped until the mean converges or until 
the noise of the entire frequency range is reached. (The frequency range between and 10 
c d _1 has a noise of 0.22 mmag.) 

b) Once the 48 frequencies were identified, the three lowest frequencies, along with their 
amplitudes and phases were fitted simultaneously by NLLS. 

c) The parameters from b) were fixed and all other frequencies, amplitudes, and phases 
were simultaneously fitted to the original light curve. 

We separate the fitting of the parameters in b) and c) to ensure convergence of the fit. 
This procedure is similar to low pass filtering of the data (i.e. removal of as much of the 
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long term trend as possible before proceeding). The parameters from (b) are necessary to 
improve the residuals. The results are presented in Table 1. Figure 2 show the DFT and the 
DFT of the comparison star and the fit parameters with bootstrapped errors (described in 
the following paragraph). The lowest frequency f3 in Table 1 has a period of 47 days which 
is longer than the data and likely to be less accurate. 

To assess the precision of our parameter-fitting method in steps (b) and (c) we use a 
bootstrap technique (see e.g. Wall & Jenkins 2003). The procedures (b) and (c) are (sepa- 
rately) repeated 10,000 times with a new light curve that is derived from points randomly 
selected from the original light curve. The new light curve has the same number of data 
points as the original with the possibility that some points may be repeated. This essentially 
introduces random gaps in the data. Figure 3 shows a typical distribution for a parame- 
ter set. The sigma of the distribution is used to estimate the 1-sigma error bars of each 
parameter. 

It is reasonable to assume that these multiple periodicities are nonradial pulsations 
simultaneously excited in the star. The analysis has revealed, in addition to the dense 
low-frequency spectrum, the presence of relatively high frequencies up to 2.8 c d _1 which 
correspond to p-mode pulsations (i.e. both p- and g-modes are excited simultaneously in 
HD 163899). To the authors' knowledge, it is the first clear detection of g-modes in a less 
luminous B supergiant, although Waelkens et al. (1998) 's list of the Hipparcos a Cyg-type 
supergiants contains some less luminous supergiants with periods of a few days indicating 
the possible presence of g-modes. 

3. MODELING THE OSCILLATIONS 

To examine the stability of post- main- sequence models corresponding to B supergiants 
and bright giants, we have computed evolution models for a mass range of 7 < M/M < 20 
with an initial chemical composition of (X,Z) = (0.7,0.02). We used a Henyey-type stellar 
evolution code with OPAL opacities (Iglesias & Rogers 1996). The stability of nonradial 
pulsations has been examined for selected models by using a finite difference method (Saio 
& Cox 1980). The stability of a selected model has also been analyzed with a Riccati shooting 
code (Gautschy et al. 1996). We have confirmed that the results of the two methods agree 
well. 

Figure 4 shows positions of models in the HR diagram which excite p-modes (triangles) 
and g-modes (inverted triangles). (The models which excite both g- and p-modes appear as 
asterisks.) If the angular frequency of pulsation is larger (smaller) than ^AnGp, the mode 
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is classified as a p-mode (g-mode), where G and p are the gravitational constant and the 
mean density of the model (see also Fig. 5). The p-mode ((3 Cephei-type) instability range 
and the usual SPB g-mode instability range for main-sequence stars of M < 9M agree with 
those obtained by Pamyatnykh (1999). We also note that the red edge of the radial-mode 
instability range of 2OM models is located at logT e g = 4.34 in perfect agreement with the 
result of Kiriakidis et al. (1993) who examined the stability of radial modes for very massive 
stars. 

In addition to the known instability regions, Figure 4 shows a wide g-mode instability 
region for the post-main-sequence massive (M > 12M ) models. The instability region 
has not been recognized before. In fact, the stability of g-modes in the post-main-sequence 
massive stars had not been examined since the emergence of OPAL opacities. [Glatzel & 
Mehren (1996) examined the stability of p-modes but not of higher-order g-modes.] In the 
next section, we focus on the the properties of unstable high-order g-modes. 

In the hotter part of the instability region p-modes are also excited simultaneously. 
The position of HD 163899 shown in Figure 4, which is roughly estimated from its spectral 
type B2 Ib/II, indicates that HD 163899 lies in the newly found instability region. In other 
words, HD 163899 belongs to a previously unrecognized type of pulsating stars, which we 
call "slowly-pulsating B-supergiants" (or SPBsg). Strictly speaking, HD 163899, showing 
both g- and p-modes, is a hybrid of SPBsg and (3 Cephei-type. 

Figure 5 shows pulsation frequencies (cycle per day) of excited modes versus effective 
temperature along the evolutionary tracks of 15M Q and 2OM . The left panel shows each 
frequency and amplitude of the periodicities of HD 163899 detected by the MOST satellite. 
The horizontal dotted lines in the middle and the right panels correspond to the MOST fre- 
quencies. The dashed lines show the frequency given by y/AjrCp/ (2tt) along the evolutionary 
track. This line is used to distinguish between p-modes and g-modes in Figure 4. Figure 6 
shows similar plots for 12M , 13M , and 14M . 

P-modes are excited in hotter models, while g-modes are excited in relatively cooler 
ones. The frequencies of the excited g-modes tend to be higher for a larger latitudinal degree 
I. In restricted T eff ranges, both p-modes and g-modes are excited simultaneously, and 
the corresponding models are appropriate for HD 163899. Since a linear stability analysis 
cannot predict the amplitudes of unstable modes, we compare in this first attempt only the 
predicted frequency ranges of excited modes with the observed ones. The observed frequency 
range is reproduced by the 2OM models around logT eff ~ 4.41 and by the 15M models 
logT c fj ~ 4.36. (These are indicated by vertical dash-dotted lines in Fig. 5.) On the other 
hand, for 14M and 13M models the frequency gap between p-modes and g-modes looks 
too large to be consistent with the observed frequency distribution. 
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Taking into account that the estimated range of logT eff for B2Ib/II is roughly between 
4.22 and 4.32 (Fig. 4), a 15M model at logT cff 4.36 looks best to reproduce the observed 
frequency range of HD 163899. Although the frequency gap is still large, it would be filled, 
at least partially, by rotational splittings. (If a higher T eff is allowed, 20M Q models at 
logTeff « 4.41 work better.) 

In our models, rotation is completely disregarded. In the presence of rotation a pulsation 
frequency measured in the observer's frame, z/(obs), shifts as 

v(obs) = z/(rot) — mfl (1) 

from the frequency in the co-rotating frame, z/(rot), where m is the azimuthal order (\m\ < I) 
and Q is the rotation frequency. 

If we assume that HD 163899 rotates at a speed of about lOOkm/s, an average for 
B2Ib/II (Abt et al. 2002), the rotation frequency would be « 0.1c d _1 . If we assume I < 3 
for observed pulsation modes, we expect that observed frequencies would be shifted by at 
most ~ ±0.3 c d _1 from z/(rot)s which we assume to be close to those obtained for Q = 0. This 
effect would nearly fill the frequency gap between p- and g-modes, and make the predicted 
frequency distribution roughly consistent with the observed one for HD 163899. 

It is not possible to perform a detailed comparison between the detected frequencies 
and theoretical predictions until accurate information on the parameters of HD 163899, such 
as the effective temperature, the luminosity, and the rotation speed is available. The most 
important finding in our modeling exercises is the fact that g-modes can be excited in post- 
main-sequence models, which suggests the presence of a new group of variable stars. We 
discuss, in the next section, the properties of excited g-modes. 

Figures 4 and 5 indicate that earlier B-supergiants should pulsate in p-modes. This 
theoretical prediction is consistent with the recent observational result of Kaufer et al. (2006) 
who detected short-period line-profile variations in HD 64760 (B0. 5 lb). 

Before we proceed to the next section, we note that there are strongly-excited very-low 
frequency modes, which are most clearly seen as a lowest-frequency branch in the frequency- 
Teg diagram for 2OM in Figure 5. The amplitude of these modes is strongly confined to the 
exterior of the bottom of the Fe-convection zone. The frequencies are so low that these modes 
are propagative even at the outer boundary. Since we use a reflective boundary condition 
at the outer boundary, we do not regard the stability of these modes to be accurate. Also, 
the property of these modes would be modified significantly in the presence of slow rotation. 
Therefore, we have not included these modes in Figure 4. We leave these modes for future 
investigations. 
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4. 



EXCITATION OF G-MODES 



We have found that g-modes are excited in massive post-main-sequence stars, whose 
central part is in radiative equilibrium. This looks contrary to the general thought that 
g-modes should be damped due to very strong dissipation expected in a radiative core where 
the Brunt-Vaisala frequency is very high. The excitation becomes possible, however, due 
to the presence of a fully-developed convective shell associated with the hydrogen-burning 
shell in a massive star. (Fully-developed convection means that convective mixing is strong 
enough to homogenize the chemical composition.) Since a g-mode pulsation is evanescent 
in the convection zone, it can be reflected at the boundary. In other words, the convective 
shell can prevent a g-mode from penetrating into the radiative core. 

Figure 7 shows runs of Lamb frequency Lf = l(l + l)c 2 s r~ 2 , and Brunt-Vaisala frequency 
N 2 , which is approximately written for an ideal gas as 



where c s is the adiabatic sound speed, r is the distance from the stellar center, g is the 
local gravity, H p is the pressure scale height, P is the pressure, and \x is the mean molecular 
weight. Frequencies of excited g-modes of I — 2 are shown by dotted horizontal lines. In this 
figure all quantities are normalized by GMR~ 3 . 

Roughly speaking, the radial wavenumber of a g-mode which pulsates with an angular 
frequency a is proportional to <y~ l \/ (N 2 — cr 2 )(Lf — a 2 ). Therefore, g-modes are radially 
propagative only if a 2 < N 2 and a 2 < L 2 , and evanescent otherwise (see e.g. Unno et al 
1989; Cox 1980, for details). Since both N 2 and Lf are very large in the radiative core, we 
expect a very large wavenumber (i.e. very short wavelength) and hence a large radiative 
damping for a g-mode in the central part. This is the reason why we do not generally expect 
g-modes to be excited in a post-main-sequence star with a radiative non-degenerate core. 

In a massive post-main-sequence star, however, a fully developed convection zone ap- 
pears at the hydrogen-burning shell, which is seen as a narrow gap around logT ~ 7.5 in 
Figure 7. This zone, where g-modes are evanescent, can reflect some g-modes and prevent 
them from penetrating into the core. This effect reduces radiative damping in the core 
significantly, and hence helps the K-mechanism in the envelope excite the mode. 

Figure 8 shows the work W, differential work dW/d(— logT), radial displacement am- 
plitude, and kinetic energy distribution (top panel) of an excited g-mode. The pulsation is 
excited if the value of W is positive at the surface. (dW/(— dlogT) > means that the 
zone drives pulsation.) We see that the amplitude and the kinetic energy suddenly becomes 
very small interior to the point at logT 7.5, where the fully convective zone is located. 
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This indicates that the reflection at the convective shell actually occurs for this mode. Such 
a reflection at the convective shell is essential for a g-mode to be excited in a post-main- 
sequence star by the K-mechanism at the Fe-bump of opacity. All post-main-sequence models 
which excite g-modes in Figure 4 have a fully-developed convection zone associated with the 
hydrogen-burning shell. On the other hand, in a stable range of 4.16 < logT e g- < 4.20 for 
13M (Fig. 6), for example, no fully-developed convective zone is present at the hydrogen- 
burning shell (i.e. the whole convectively-unstable layers are semiconvective only). 

Only a few selected g-modes are reflected at the convection zone and excited by the Fe- 
bump in an appropriate frequency range. Figure 9 compares the kinetic energy distribution 
of an excited I = 2 g-mode (thick solid line) to that of a damped mode (thin solid line) as a 
function of the fractional radius (x = r/R). (In this diagram the right hand side is the stellar 
surface.) The pulsation frequencies of the two modes are very close; 0.324 c d -1 and 0.326 
c d _1 . These eigenfunctions were calculated by the Riccati shooting method (Gautschy et 
al. 1996), completely resolving the spatial oscillations with roughly 16,500 integration steps. 
It is obvious that the kinetic energy in the core of the excited mode is significantly reduced 
and hence the radiative damping is quenched, while the kinetic energy of the damped mode 
is strongly confined to the central part. 

Most of the g-modes in the frequency range where excited modes reside are confined 
to the radiative core and are damped due to strong radiative damping. Only those selected 
g-modes whose spatial oscillations have the right phase at the boundary of the chemically 
homogeneous shell convection zone are reflected and confined to the envelope. Therefore, 
which g-modes are excited depends on details of the structure around the convective shell. In 
the model shown in Figure 7 the layer in 7.23 < logT < 7.44 is semi-convective (i.e. convec- 
tive zone with an inhomogeneous chemical composition). We have treated semi-convective 
mixing based on the model of Spruit (1992). Everyone admits that many uncertainties lie 
in the theory of convective mixing. Since the Brunt- Vaisala frequency depends sensitively 
on the distribution of chemical elements [Eq.(2)], frequency spacings of the excited g-modes 
depend on the treatment of convective mixing. This indicates that comparing observed fre- 
quencies with the predicted ones from models can constrain the theory of convective mixing. 
We also note that the lower luminosity boundary of the g-mode instability region, which is 
the boundary of the presence of a fully developed shell convection, depends on the efficiency 
of convective mixing. Thus, g-modes in SPBsgs have great potential for asteroseismology. 

As discussed above, the presence of a convective shell is crucial to excite g-modes in the 
post-main-sequence models. The appearance of a convective shell in the post-main-sequence 
phase of massive stars has been known for more than 30 years (e.g. Barbara, et al. 1972). 
It has also been known for a long time that g-modes can be reflected at the boundary of 
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a convection zone, and that the reflection quenches a strong radiative dissipation in the 
core. Nevertheless, we can find no published evidence of stability analyses of g-modes in 
massive post-main-sequence stars since the emergence of OPAL opacities. The discovery of 
light variations in HD163899 by the MOST satellite, which prompted us to perform stability 
analyses for such modes, played a crucial role in our eventually finding the excited high- 
order g-modes. On the HR diagram, the region with unstable g-modes seems to extend to 
the region of a Cyg variables suggesting that g-mode pulsations might play an important 
role in understanding these luminous variable stars. The spatio-temporal evolution of the 
H-shell convection zone in very massive stars during their post-main-sequence evolution can 
be judged from figures in e.g. Lamb et al. (1976) or Maeder (1981). 



5. CONCLUSION 

The MOST satellite discovered that the B type supergiant HD 163899 pulsates in g- 
modes as well as in p-modes, showing, for the first time, SPB-type g-mode pulsations in a 
less luminous B supergiant. We have performed pulsation stability analyses for evolutionary 
models in the mass ranges of 7 < M/M Q < 20, and found that g-modes are excited by the 
Fe-bump of opacity in post-main-sequence stars in a wide range of effective temperatures if 
the stellar mass is larger than ~ 12M . Our results indicate the existence of a previously 
unrecognized group of variable stars; slowly pulsating B supergiants (SPBsg) distinct from 
a Cyg variables. 

The excitation of g-modes in B supergiants becomes possible by a reflection of g-modes 
at the fully-developed convection zone associated with the hydrogen-burning shell. This 
reduces significantly the radiative damping in the radiative core of the star. Since the 
frequencies of excited g-modes depend on the structure around the convective shell, SPBsg 
have a great potential for asteroseismology. 

HD 163899 is the first member of the SPBsg family. Strictly speaking, HD 163899, 
showing both g- and p-modes, is a hybrid of SPBsg and (3 Cephei-type. The B supergiant 
l CMa (B3 lb /II) could be another SPBsg member. Balona & Engelbrecht (1985) obtained 
a 0.717 c d _1 periodicity for i CMa. Since the light curve looks irregular, it is likely multi- 
periodic. Another possible candidate could be HD 98410 (B2/B3 Ib/II) which is listed in 
Stankov & Handler (2005) as one of the "poor and rejected (3 Cephei candidates" with 
a Hipparcos-deduced period of 1.453 days. We also note that Waelkens et al. (1998)'s 
list for a Cyg-type supergiants includes some less luminous B supergiants or bright giants 
such as HD54764 (Bllb/II), HD141318 (B2II), and HD168183 (Bllb/II). These stars could 
also be SPBsg candidates. Further observations of these stars are needed to confirm their 
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membership. Also, for a detailed asteroseismological study of HD 163899, spectroscopic 
observations are important to properly identify the pulsation modes that are responsible for 
the photometrically observed frequencies. 
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Table 1. HD 163899 periodicities from MOST photometry 
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Table 1 — Continued 
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a Phases are referenced to the first observation in the data 
set. 
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Fig. 1 — a. The 2005 MOST 37 day light curve of the guide star HD 163899 and the 
simultaneously observed guide star HD 164388 whose light curve is shifted by 30 mmag for 
better visibility; b. a 2 day portion of the HD 163899 light curve with the data binned 
every 2 min. The 2a errors are ~ 1.3 mmag. Because MOST was operating outside the 
continuous viewing zone the duty cycle was ~50% accounting for the gaps in the data. The 
solid line in the bottom panel is the fit of the 48 frequencies detected to the data. A short 
term variability of ~ 0.05d seen in both stars arose from the orbital phase of the MOST 
satellite, and was not included in our frequency search. 
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Fig. 2. — Amplitude spectra from discrete Fourier transforms for the light curves of HD 
163899 (solid line) and the comparison star HD 164388 (dashed line), where the sign of the 
amplitude of the comparison star is inverted for better visibility. Dots and triangles with 
error bars indicate the frequency/amplitude obtained for HD 163899. Peaks at 1 and 2 c d _1 
conspicuous in the comparison star are artifacts typical of the MOST photometry. 
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Fig. 3. — The bootstrap distributions of parameters for the frequency fi. They represent the 
number of the cases at each parameter bin resulted from 10,000 light curves produced by 
inserting random gaps into the original. At the top of each diagram the mean value \i and 
the standard deviation a of the distribution are indicated. Each diagram from the bottom 
to the top shows the distribution for the frequency v in c d" 1 , the amplitude in mmag, and 
the phase in radian. 
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Fig. 4. — Massive star models which excite p-modes (triangles) and g-modes (inverted trian- 
gles) are shown along the evolutionary track for each mass. The models which excite both 
g- and p-modes appear as asterisks. A large square with error bars indicates the position of 
HD 163899 estimated from the spectral type B2 lb /ll. 
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Fig. 5. — The frequencies of excited pulsations of < I < 3 in evolutionary models of 
15M Q and 2OM against the effective temperature. The left panel shows frequencies and 
amplitudes of pulsations obtained by MOST for HD 163899. The horizontal dotted lines in 
the middle and the right panels correspond to the MOST frequencies. Long dashed lines 
in the middle and the right panels indicate the relation of a — y/AnGp, which is used to 
discriminate between p- and g-modes in Fig. 4, where a is the angular frequency of pulsation, 
G is the gravitational constant, and p is the mean density of the star. A vertical dash-dotted 
line in each of the middle and the right panels indicates a value of log T e s around which the 
frequencies excited in the models are roughly consistent with the observed frequencies. 
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Fig. 6. — The frequencies of excited pulsations of < / < 3 in evolutionary models of 12M Q , 
13M and 14M against the effective temperature. The dashed lines and dotted lines have 
the same meanings as those in Fig. 5. 
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Fig. 7. — Runs of the square of the Brunt- Vaisala frequency, N 2 , and the square of Lamb 
frequency Lf = 1(1 + l)c 2 s /r 2 for I = 1,2,3. Horizontal dotted lines indicate the angular 
frequencies (a) of g-modes of I — 2 excited in this model. All the quantities in this figure 
are normalized by GMR~ 3 . 
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Fig. 8. — The bottom panel shows radial displacement (solid line), work W (dotted line), 
and differential work dW/(— dlogT) (dots), and opacity k (dash-dotted line) as a function of 
log T. The top panel shows the distribution of the kinetic energy of the pulsation. Radiative 
dissipation in the core is quenched because the pulsation is reflected at the convective shell 
at logT PS 7.5. The opacity peak at logT ps 5.3 arises from Fe and that at logT ps 4.6 from 
Hell ionization. 
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Fig. 9. — Kinetic energy distribution of an excited 0.3241 c d^ 1 mode (thick solid line) and 
a damped 0.3256 c d _1 mode (thin solid line) of / = 2 as a function of the fractional distance 
from the center (x = r/R). (The thin line is the upper one.) Dot-dashed lines show the 
suitably scaled spatial variation of the mean molecular weight. The kinetic energy of the 
excited mode is very small in the core due to the reflection in the convective shell, which 
reduces radiative damping in the core and hence makes the mode unstable. 



